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The possibility of extending the Liouville Conformal Field Theory from values of the central charge 
c > 25 to c < 1 has been debated for many years in condensed matter physics as well as in string 
theory. It was only recently proven that such an extension—involving a real spectrum of critical 
exponents as well as an analytic continuation of the DOZZ formula for three-point couplings—does 
give rise to a consistent theory. We show in this Letter that this theory can be interpreted in terms 
of microscopic loop models. We introduce in particular a family of geometrical operators, and, 
using an efficient algorithm to compute three-point functions from the lattice, we show that their 
operator algebra corresponds exactly to that of vertex operators Vk in c < 1 Liouville. We interpret 
geometrically the limit d —>■ 0 of Va and explain why it is not the identity operator (despite having 
conformal weight A = 0). 

PACS numbers: 05.70.Ln, 72.15.Qm, 74.40.Gh 


Quantum Liouville theory has played a fundamental 
role in string theory and quantum gravity since its in¬ 
troduction by Polyakov [T]. It also provides a crucial, 
exactly solvable, example of non-rational conformal field 
theory (CFT) [2]. These theories are difhcult to study, 
and ripe with unusual features, such as continuous spec¬ 
tra of critical exponents. They are however believed to 
play an important role in many applications, the best 
known being the transition between plateaux in the in¬ 
teger quantum Hall effect laa. 

Recent work on Liouville theory has focussed on con¬ 
tinuations of the parameters to new domains of values, 
where new applications may be found, while the exten¬ 
sion of the standard CFT solution presents considerable 
challenges. We will focus in this Letter on the case of 
central charges c < 1, with the Liouville action 


A = 



-b iQUcj) + 47r/xe-2*^'^ 


( 1 ) 


where 5“^ is the metric and TZ the Ricci scalar of the 
underlying 2D space. The coupling constant b is real, 
and Q = {b~^ — b). The corresponding central charge 
is c = 1 — 6(3^, and the conformal weights of the vertex 
operators Va = 6^““^ are 

A = A = d(d - Q). (2) 


We note that this CFT is often called “time-like Liou¬ 
ville” , since, under a redefinition (j) = i(j), the dynamics 
for (j) looks like ordinary Liouville except that the kinetic 
term has negative sign, just like the time coordinate in 
Minkowski metric. However, since a certain domain of 
values of a is usually implied in time-like theories, we 
instead refer to 0 as “c < 1 Liouville”. 


The extension of Liouville theory to the case c < 1 was 
proposed in several contexts. A chief motivation was the 
hope to compute correlation functions in non-minimal 
models of statistical mechanics [S] — for instance, in the 
percolation problem with c = 0 [6]— which have remained 
largely unknown up to this day. An essential step in this 
extension was the discovery of a formula for the three- 
point function of vertex operators, generalizing the so- 
called DOZZ formula for ordinary Liouville dll], which 
was the starting point of the construction of the corre¬ 
sponding non-rational CFT. This formula reads [Hilo] 


A,. . AgT(6-Q + di23)nT(S + d5) 

6(01,0:2,0:3) = ! = , 

VnLiTr(^ + 2aOT(S-Q + 2d,) 

(3) 

where the product J([ makes {ijk) run over the three 
cyclic permutations of (123), = ai -\- aj — and 

di23 = di -b 02 + 0:3. The normalization condition 
C'(d,d,0) = 1 defines Ag. The function C'(di,d2,d3) 
is totally symmetric in its three arguments, and is also 
invariant under di —>■ Q — di for any i. The function T 
is given by (setting q = b -\- b~^) 


lnT(a;) 




sinh^ (f -x)| 
sinh ^ sinh 

^ 2b 


(4) 

for 0 < Re(x) < q. Outside that range, T is extended us¬ 
ing functional relations; see the Supplementary Material 
(SM). C{ai,a 2 , d3) is believed to encode the three-point 
function of vertex operators, whose general form is die- 
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tated by conformal invariance 

{V&Ari)V&^{r2)Va^{r:i)) = C(di, 02, , ( 5 ) 

where Vij = jr^ — r^j. The two-point function of opera¬ 
tors with identical values of a is normalized to have unit 
residue. 

The consistency of the c < 1 theory with a = ^ +p and 

p G M—viz. A = p'^ — ^ by ([^ —was recently demon¬ 
strated in m, including by extensive numerical checks 
of crossing symmetry for the corresponding four-point 
functions. Meanwhile, the relevance of Liouville theory 
to conformal models of fluctuating loops was pointed out 
in [laiii!. These works, based on the so-called geomet¬ 
rical Coulomb gas (CG) construction [ 14 ], were however 
limited to two-point functions. But following the sug¬ 
gestion in |S], an interesting proposal was made that 
the probability for three points to belong to the same 
Fortuin-Kasteleyn (FK) cluster in the critical Q-state 
Potts model was simply related to the three-point cou¬ 
pling (§ for a particular value of the charges cti = a{Q). 
This was confirmed by numerical simulations m for real 
Q G [ 1 , 4 ], up to a factor ^/ 2 . In | 5 |, this factor was 
explained by relating the correlation functions of spin 
operators to those of kink operators, through Kramers- 
Wannier duality. Due to the permutation symmetry of 
the Potts model, the kinks exhibit a two-channel struc¬ 
ture in their OPEs, which produces a factor in any 
measurement of the structure constant. 

In this Letter we provide a complete statistical physics 
interpretation—and extensive numerical checks—of the 
three-point coupling ([^, for continuous values of c and 
the three independently varying “electric” charges oti. 
Our results apply both to the loop model underlying the 
Q-state Potts model, and to the dense and dilute phases 
of the 0 (n) loop model [TB]. In the case of the Potts 
model, we confirm the factor found in , and we re¬ 
cover a two-channel structure at the level of the transfer 
matrix Hilbert space. 

The context of this physical interpretation is the Con¬ 
formal Loop Ensemble (CLE) [ 17 ] . familiar in the study 
of spin models, CG mappings UlllH], and more recently 
the Schramm-Loewner Evolution (SLE) [TSj. The most 
physical lattice discretization of this ensemble consists 
m in drawing self- and mutually avoiding loops on the 
hexagonal lattice, with a fugacity per loop equal to a real 
number n G [— 2 , 2 ] and a fugacity per monomer (3 which 
is / 3 e = (2 + (2-n)i/2)i/2 [ig in the dilute phase, and 
/ 3 c < /3 < 00 in the dense phase; both phases are critical. 

To link with the CG approach, we parametrize n = 
—2cos7rg, with coupling g G ( 0 , 2 ], where g G ( 0 , 1 ] (resp. 
g G [ 1 , 2 ]) describes the dense (resp. dilute) phase. Set¬ 
ting eo = 1 — 5 the central charge is c = 1 — beg//? [H] . 
This matches c < 1 Liouville, provided b = and the 
link to CLEk is k = 4 ,/g. 


The question of which CFT describes the CLE com¬ 
pletely remains open to this day. This CET should in¬ 
volve a free boson, which one can interpret as the long- 
wavelength limit of a solid-on-solid model (SOS) dual 
to the loops. More specifically | 20 |, the weight n per 
loop must be understood as arising from oriented loops, 
which get an extra (complex) weight for every left and 
right turn. On the hexagonal lattice, this weight is 
exp(±i 7 reo/ 6 ) so that, after summing over both orienta¬ 
tions, and using that the number of left minus the number 
of right turns of a closed loop on the hexagonal lattice is 
±6, produces the correct fugacity n = 2 cos ttco per loop. 
The oriented loops define the SOS model by duality, 
heights on neighboring faces differing by ±Ah depend¬ 
ing on the arrow that separates them. Finally, the fact 
that the weight depends on the number of turns can also 
be construed [Hdl] to explain the second, curvature- 
dependent term in Q. 

The key result of this CG analysis is the understanding 
of two-point functions. Specifically, the partition function 
of the loop model with two special points fq, rq—such 
that loops separating rq from fq get the fugacity ni = 
2 cos7rei (with ei arbitrary) instead of n —decays, for 
large distances, as with 

2 2 

A(ei) = ^i^ = di(ai-g). (6) 

In the second equation, we have formally matched the 
well-known CG result m with the c < 1 Liouville for¬ 
mula ([^ , suggesting that the (non-local) observable mod¬ 
ifying the weight of the loops is related with a vertex 
operator of charge 

ai = ?+P, withp=^. ( 7 ) 

The symmetry di —> Q — di then amounts to ei -G — ei. 

The geometrical CG approach is equivalent, in GFT 
parlance, to a free boson with a charge at infinity [22] . 
Up to this day, no general result has been available for the 
three- and higher-point correlation functions in the loop 
model. While part of the difficulty lies in determining 
the correct definition of these correlations in geometrical 
terms, there are also deep conceptual issues to be over¬ 
come. For instance, many non-trivial correlations seem 
to exist, which the screening construction in |22j would 
erroneously set to zero for reasons of charge neutrality. 

Inspired by the discussion in um and the observa¬ 
tion in [adg we have investigated the possible meaning 
of the three-point function ([^ within the loop model. 
We have found overwhelming evidence that it admits a 
geometrical interpretation, which is as follows. 

Consider three points rq, fq, fq in the plane, and imag¬ 
ine we run a cut C12 from fq to fq, and another cut 
C23 from 7*2 to fq. We now define a modified partition 
function ^ni,n2,n3(fi, f^2, f^a) of the loop model by giving 



?! 



FIG. 1: Loop weights in the three-point function, (a) Generic 
case. Black loops have the bulk weight n, while red, blue and 
green loops have weight m, n 2 and na respectively. The figure 
shows all topologies simultaneously; in reality loops cannot 
intersect, (b) The same, but with point ra sent to inhnity. 


four different weights to the different topological classes 
of loops (see Fig. If Ni 2 = 0,1 (resp. N 2 ^) denotes 
the number of times modulo 2 that a given loop intersects 
Ci 2 (resp. (723), then its weight is 

rii , with i = 2 N 23 fVi 2 G {0,1, 2, 3} . (8) 


We have here set no = n. In other words, a loop sepa¬ 
rating fi from the other two points gets weight rii (with 
i = 1,2,3), while a loop surrounding none or all three 
points gets the bulk loop weight Uq. 

Notice that since a loop can be turned around the 
“point at infinity” on the Riemann sphere, we cannot 
distinguish a loop surrounding a subset of points from 
a loop surrounding its complement. Therefore iV-point 
functions allow for 2^“^ distinct weights. To weigh dif¬ 
ferently all 2^ ways of surrounding subsets of N points, 
we need to consider an (N + l)-point function with tn+i 
sent to infinity (in particular no loop can surround fV+i). 
This is shown in Fig. & for N = 2. 

Parametrizing the weights rii = 2 cos Tre^ with Cj G 
[—1,1], and using 0 provides a set of three charges 
with * = 1, 2,3 (see SM for higher values of e^). The key 
idea is that the partition function with weights (|^ is pro¬ 
portional to the three-point function of the vertex oper¬ 
ators Vai in c < 1 Liouville. To make this statement def¬ 
inite we need to impose the correct normalization of the 
partition function. For conciseness we render implicit the 
insertions r) and abbreviate Zni,nj,nk{Ti,r 2 ,r^) = Ziji^. 
We then have our main result that .^123 is proportional 
to ([^, and more precisely 

C{ai, Q! 2 , ds) = 


^123 






oil 


^202 


^330 


000 " 


■^lOl^^llO ^'220^022 Z033Z303 ^ ^ 

The normalization on the left-hand side corresponds to 
setting C(d, d, 0) = 1. 


Dense 0(n) model Dilute 0(n) model 



FIG. 2: (7(d, a, a) as a function of ni = n 2 = ns in the dense 
and dilute 0(n) model with n = 1. 


In order to check these formulas, we have devised a 
method to determine three-point couplings numerically 
using transfer matrices (see SM for more details). It 
turns out more convenient to study models defined on the 
square lattice, rather than the honeycomb 0 {n) model 
[l 6 ] discussed above. The square-lattice 0(n) model 
P5] has again dilute and dense phases (representing the 
same universality classes). Moreover, we studied the in¬ 
state Potts model via its equivalent completely packed 
0 (n = ^/Q) model on the square lattice [3D], which pro¬ 
duces only the dense universality class. 

In our numerical scheme, the axially oriented square 
lattice is wrapped on a cylinder of circumference L. We 
split the cylinder into two halves, each consisting of 
M ^ L rows, and place ri,r 2 ,r 3 at the bottom, middle 
and top respectively, all at identical horizontal positions. 
The boundary conditions at the bottom and top are cho¬ 
sen such that neither fi nor rs can be surrounded by a 
loop. Then ([^ amounts to giving weight ni (resp. 713 ) to 
non-contractible loops below (resp. above) r 2 , and weight 
712 to contractible loops that surround 7 ( 2 - All other loops 
get weight n. We then determine numerically the corre¬ 
sponding partition functions by acting with the transfer 
matrix, and form the ratio corresponding to Q after the 
conformal map from the plane to the cylinder: 


C{ai, 0.2, 03) 


rycyl / vcyl rT’cyl 
^123 / -^202 -^000 

ycyl \ ycy\ ycy\ ’ 

^220 V ^101^303 


( 10 ) 


where Z^jj, is the analog of on the cylinder, with 
the points fq and fq sent to the two extremities of the 
cylinder. We finally compare the values of this ratio with 
the Liouville result ([^, evaluated as described in SM. We 
have found excellent agreement for a wide range of values 
of the four loop weights. We shall now illustrate this by 
reporting in detail on two cases which are particularly 
noteworthy. 
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The first case corresponds to ni = 712 = na, for which 
Fig. shows results both for the dense (c = 0 , perco¬ 
lation) and dilute (c = |) phases of the 0 (n) model 
with n = 1 . The agreement with (|^ is excellent for 
the whole range of rii, and can be further improved by 
finite-size scaling (FSS) extrapolation (see SM). The di¬ 
vergence C —>■ —00 as rii —>■ — l+j due to a pole in the 
theoretical formula, is beautifully reproduced by the nu¬ 
merics. Note that to reach values near = — 1 requires 
analytical continuation of T beyond the domain of valid¬ 
ity of the integral representation Q; see SM for details. 
The formula (|^ works perfectly well for > 2 as well, 
corresponding to imaginary values of the charges a. We 
conjecture that ^ applies indeed for all Ui S C, corre¬ 
sponding in general to complex values of u,. For n S R, 
we need — 2 < n < 2 to ensure a critical theory; nothing 
is known about complex values of n. 

We stress that for generic (ni,n2,n3), the partition 
function „3 cannot be encoded in the local vertex 

model equivalent m to the 0 (n) model. Indeed, in the 
former, Zn-^,n2,n3 can only be obtained by introducing a 
twist factor e^^®^ (resp. associated to the arrow 

flux through the cut C12 (resp. €23)^ and this forces 712 
to take one of the four values —2cos7r(g ± ei ± 63). In 
the CG approach | 52 ], this corresponds to a three-point 
function satisfying the neutrality condition with one type 
of screening charges. Hence, going beyond this case is 
only possible in the non-local loop model. 

Repeating the numerics for the loop model associated 
to the Potts model, we still find excellent agreement for 
Ui > 0 . For Tii < 0 we have very large FSS effects, but 
going to large sizes {L = 16 ) the extrapolation still gives 
a decent agreement with the theoretical formula. Right 
at Tij = 0, and only there, one can check with great ac¬ 
curacy that the numerical data converge to \/ 2 C (and 
not to C as usual). This is precisely the case 0 ITS] , 
where ([^ is interpreted as the probability that rq, 7^2, 
belong to the same FK cluster. Actually, one can show 
[211 [2S] that the space of states on which the transfer ma¬ 
trix is acting splits, when = 0 and only then, into two 
isomorphic subspaces, which are technically irreducible 
representations of the lattice (periodic Temperley-Lieb) 
algebra underlying the dynamics of the model; see SM for 
more details. Numerical methods and their associated 
normalizations measure the three-point constant within 
one subspace only (selected by the boundary conditions 
imposed at the ends of the cylinder), while, by analytic- 
ity, the Liouville result holds for the whole space of states. 
This leads, by easy considerations on ( 10 1, to a measured 
three-point constant which is too large by a factor for 

Hi = 0 . 


The second case we consider in some detail here is 
when one of the charges vanishes, say 0:2 = 0 . The vertex 
operator then has = 0 . In ordinary CFT this would 
imply 14,2(7^2) = /, the identity operator, so that (§ 
reduces to a two-point function, equal to zero by confor¬ 


mal invariance, unless di =0:3. However, ([^ does not 
exhibit this feature at all, and the function C remains 
highly non-trivial even when one of the arguments van¬ 
ishes |2S]- We still have C'(d, 0 , d) = 1 , compatible with 
the normalization of two-point functions, but in general 
C'(di,0,d3) ^ 0 . This will be illustrated numerically 
below. 

To discuss the corresponding geometrical meaning, we 
consider the loop model with 712 = n. According to 
Fig.^ the weight of a loop encircling rq (or rq) depends 
on whether it also encircles fq- Thus, fq is not invisible 
at all, and Vq2(T2) is in fact a marking operator, distinct 
from /, even though Aq,^ = 0 . The three-point function 
([^ is then 


■^103 

Z 


C{ai, 0, ds) 



/ X Aa, -l-Aa, 


where we have, for the time being, considered an unnor¬ 
malized quantity, where the (unknown) lattice cutoff a 
appears explicitly. This depends on Tq (despite A^^ = 0 ), 
implying a non-zero derivative wrt fq. Therefore L_i 
acts non-trivially on and this field is non-degenerate. 
We now send fq —)■ cx). The fraction of loops encircling 
both rq and at least one of the other two points becomes 
negligible, so we get rid of the marked point. The first 
scale factor above disappears, and we can then get rid of 
the second one by normalizing with the two-point func¬ 
tion. Hence 


Zcy\ 

ni,n3 




ni 


C , 0,0^3) 5 


(11) 


where now denotes the partition function where 

loops encircling fq (resp. rq) get weight 711 (resp. 713), 
while those encircling both points get weight n. 

Alternatively, we can keep 712 generic (so that Aq^ ^ 
0 ), send rq —>■ 00 as before, but choose the microscopic 
boundary conditions so that no loop can encircle 7 q- 
Loops encircling both rq an d ry then get the non-trivial 
weight 772 , and the ratio ( [ll| produces C{ai,a2,c(3). 
This situation is depicted in Fig. [^. 

We can also simplify © even further by noticing that 
C remains non-trivial even when two of the charges at 
vanish. We have then 


Zcyl 
ni ,n 


vcyl 7cyl 


= ( 7 (ai, 0 , 0 ). 


( 12 ) 


Here, is the partition function where loops encir¬ 

cling 7 q but not 7 q get a weight ni. 

We have checked these interpretations of C(di,0,d3) 
and C'(q;i, 0 , 0 ) numerically. The latter case is shown in 
Fig. [31 For the Potts model one obtains similar results, 
except that now C/v^ is observed when 77i = 0, as can 










5 


Dense 0(n) model Dilute 0(n) model 



FIG. 3: < 7 ( 0 : 1 ,0,0) and (5(0 ,q:i, 0) as functions of ni in the 
dense and dilute 0(n) model with n = 1. The two positions 
of the vertex operator Vax —at one extremity or in the middle 
of the cylinder—give different microscopic results, but their 
L —>■ 00 limits agree with the same analytical formula. 


again be seen from (10) and the splitting of the state 
space. 

Summarizing, the loop model with four different 
weights (see Fig. gives a consistent interpretation of 
the three-point function in c < 1 Liouville. It is also 
consistent with interpreting C as an OPE coefficient. 
This can be seen geometrically by considering the limit 
r ’23 —>■ 0. The amount of loops with weight 712 or is 
then negligible, since they would have to be “pinched” 
between r 2 and ffi. Instead, loops get weight ni when¬ 
ever they encircle ri or r 2 , and n if they encircle both. In 
other words, we recover—up to a numerical factor, and 
after subtracting the divergence at r 2 = —the two- 
point function of operators with charge di, which ex¬ 
plains in what sense a continuum of di charges appears 
in the fusion of charges d 2 and da. 

The loop model also gives natural explanations to a 
well-known paradox encountered in trying to make c < 1 
Liouville into a consistent CFT. We see here that the 
vertex operator with d = 0 must be interpreted not as 
the identity but as a ‘marking’ operator, a feature very 
similar to what happens in the SLE construction. For in¬ 
stance, in the CFT proof m of Schramm’s left-passage 
formula |28j a point-marking (or “indicator” [ 211 ) op¬ 
erator of zero conformal weight—but distinct from the 
identity operator—is used to select the correct conformal 
block in the corresponding correlation function. Similar 
arguments can be made within boundary CFT |29] . 

We emphasize (more details in SM) that nothing spe¬ 
cial seems to appear as c —>■ 1. In particular, the loop 
model does not produce the singularities in the three- 
point function obtained in [30] . This is, from our point of 
view, quite expected, since the model in [30] is obtained 
as the m —>■ 00 limit of the Am RSOS models, whose 
relationship with the loop model is not straightforward. 


and involves a complicated sum over sectors. 

We hope to report later on the study of iV = 4 point 
functions, and to explore the geometrical definitions nec¬ 
essary when even more different loop weights are intro¬ 
duced. Boundary extensions of this work also offer tan¬ 
talizing possibilities. 

Similar three-point functions can be defined involving, 
instead of loops with modified weights, the so-called wa¬ 
termelon or fuseau m operators. In CG parlance, the 
electric operators Eq. are then replaced by magnetic op¬ 
erators Oe- inserting £i defect lines. We have determined 
numerically the corresponding C, the simplest of which 
(all ii = 2) determining the probability that three points 
lie on the same loop. However, the results do not match 
formula ([^ when the corresponding Ai. are entered, and 
seem to correspond to another continuation of the Li- 
ouiville results than the one considered here. This is but 
one of the many aspects (such as the physical meaning 
of the pole in the three point function) of this problem 
that deserve further study. We hope to report on all of 
this soon (some related results have appeared in [32]). 
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SUPPLEMENTARY MATERIAL 
The function T(x) 

We first discuss the special function T (x) appearing in 
the expression ^ for the structure constants C. Intro¬ 
ducing q = b + b~^ (not to be confused with Q), it is 
related to the double gamma function r 2 (a:|a;i,a; 2 ) via 

T(x)= - r, T-^ix)=T,ix\b,b-^). (13) 

rg(a;)rg(g-a;) 

It obeys the reflection property 

T{x) = T{q-x) (14) 

together with the functional relations 
T(g/2) = 1, 

T{x + b)=-f{bx)b^-^'^^ T{x), (15) 

T(a; -f b-^) = j{xb-^) T(a;), 

where j{x) = r(a:)/r(l — x) and F is the Euler gamma 
function. 

The integral definition Q of T(x) converges only if 
\q — 2 Re(a:)| < q, i.e., 0 < Re(a:) < q. Outside this 
range, the function has to be extended using the foregoing 
formulae. 

Note that T(a;) is entire analytic with zeros at 
X = —mb~^ — nb , 

(16) 

x = {m+l)b ^-I-(n-I-1)6, 
for any m,n G Z-^. 


Special cases of C(ai, 0 : 2 , 03 ) 


In the main text we have discussed in details two spe¬ 


cial cases of our principal result (10). In the first case 


(all Hi equal), the structure constant is 


Consider for instance the structure constant ( [T7| for 
a bulk value n = 1, as plotted in Fig. As the weight 
of the special loops Ui is lowered from = 2 , a zero 
of the second T function in the numerator is first en¬ 
countered for 2 — = 0 or e = corresponding to 

i 2 d 


Ui = 2cosl^ = 2cos^ = —0.34. Beyond this value, 
the naive evaluation of the integral gives InT = — 00 , so 
the naive structure constant is zero. This is overcome by 
using the second relation in (15) and thus replacing 




?3eb 


r 1 

f i-b-^ 

1 3eb-2^ 

1 ' 

K 2 

+ 2 ) 

ri 

(i±|^ 

3eS-2 \ 
2 ) 

• ^ 

, 


6-f 36-^ 

3e6 ^ ^ 


2 

2 


(19) 


Another remark concerns the regime rii > 2. Although 
the integral (|^ is convergent here, its naive evaluation us¬ 
ing standard software such as Mathematica or Maxima 
poses problems, because the integrand decreases fast with 
t. The practical solution is to constrain the numerical in¬ 
tegration to a suitably chosen finite interval, t G [ 0 , tmax], 
outside which the integrand is exponentially small. 


Limit of central charge c = 1 

The formulae simplify considerably in the limit c = 1, 
and allow a comparison with the construction of [30]. 
The electric charge is then defined by e = 2d, and the 
conformal weights of the vertex operator read: 

A = A = d^ = ^ . (20) 

4 

In the regime when 0 < < 1 and 61 - 1 - 62 - 1-63 < 2, 

the arguments of T in ([^ all lie in the interval ]0, q = 2[, 
where the integral representation ([^ is valid, and after a 
change of variable /3 = e“‘ in 0 , one obtains 

C'(di,d2,d3) = exp[(3(ei,e2,e3)], (21) 


where 


C'(d, d, d) oc 


T3/2 (5-1(1 _e))T3/2(S-i(l + e)) ' = f 

Jo 


d/3 


o(-ln/3)(l-/3)2 


2 +^ (/3“ + r 2 


In the second case (only one rii ^ n) we have 


2 2 
T 


f 3b-b~^ I eb~^ A 

V 2 + 2 ; 


TI/2 (6 + e6-i) Ti/2 (5-e6-i) 


Because of the finite domain of convergence of the inte¬ 
gral, analytic continuation is necessary to obtain results 
for large ranges of the lattice parameters. 


£=±1 


+/3n - E p(€^x+eyy+€;^z)/2 


( 22 ) 


The structure constants for other values of the electric 
charges can be related to the above regime through the 


■ (18) functional relations (15). Note also that for 6=1, the 


double gamma function in (13) is related to the Barnes 
G-function by 


r 2 (a:|l, 1 ) = const x 


(27r 


l ®/2 


G{x) 


(23) 





























0(2) model Detail 



FIG. 4: C{a., a, a) as a function of n\ — n 2 = ns in the dense 
0 (n) model with n = 2, compared with the c = 1 limit of §. 


overlap vanishes when the radius of the cylinder L —> oo 
as a power law of L, a manifestation of the well-known 
Anderson orthogonality catastrophe [33]. However, this 
power law is factored out by the normalization (see be¬ 
low), which extracts directly the corresponding ampli¬ 
tude. Of course, this amplitude is trivial: one can check 
that C(ai, di -|- da, da) = 1, since Q(x, x + z,z) = 0, as 
follows from the simple identity (^ = — In /3) 


We further emphasize that the C function obtained in 
the loop model does not exhibit any of the singularities 
observed in the limit of minimal models studied in [30]: 
everything behaves as if the prefactor P{x,y,z) in this 
reference were simply absent. 


The case c = 1 is particularly interesting since the loop 
model then maps, in the standard Coulomb gas formal¬ 
ism [Hdl], to a free boson with no electric charge at 
infinity. This is because loops in the bulk get a weight 
n = 2, which can be obtained simply by summing over 
two possible orientations. The standard phenomenol¬ 
ogy m then does not suggest the presence of “floating 
charges”. Nonetheless, the Liouville description is again 
fully confirmed by the numerical study of the loop three- 
point functions. 

We illustrate this by considering the first case of three- 
point couplings for identical loop fugacities, C'(d,d,d). 
As seen in Fig. [^ at the point ni = n 2 = ns = — 1, 
corresponding to ei = 62 = 63 = 2/3 where (22 1 diverges, 
the structure constant observed numerically in the loop 
model behaves smoothly, as predicted by the Liouville 
description. This contrasts with the results of [30] , where 
the unitarity constraints in minimal models are reflected 
as a non-analytic behaviour of C in the limit c —)■ 1, 
which would produce (7 = 0 all along the range — 2 < 
ni = 712 = u-s < — 1. This particular example shows that 
the c —?► 1 limit of the Liouville action 0 and that of 
minimal models lead to different theories, and that the 
loop model clearly relates to the Liouville theory. 


Orthogonality catastrophe 

A potential source of confusion must be dispelled at 
this stage. Indeed, the way we determine numerically the 
three-point function on the cylinder is reminiscent of the 
calculation of overlaps of ground states for twisted spin 
chains. While one has to be careful with this correspon¬ 
dence in general, it becomes exact in the case n = 2 and 
0(2 = di -|- 0:3, where the partition function used to de¬ 
termine C is proportional to the overlap of XXX ground 
states with two different twists. When di ^ 0 ( 3 , this 


Loop models on the square lattice 


The loop model on the square lattice has nine possible 
configurations at each vertex: 



Each vertex has a Boltzmann weight pi as indicated, 
and there is a weight n for each closed loop. In the com¬ 
putation of three-point functions, some of the loops will 
get modified weights rii, 77,2 or 773 (see below). 

The square-lattice 0 ( 77 ) model referred to in the main 
text corresponds to a set of integrable weights [33] which 
read, at the isotropic point, 

Pi = 1 -t- sin A -I- sin(3A) — sin(5A), 

P 2 = P 3 = P4 = P 5 =‘2' sin(2A) sin((6A -I- 7r)/4), 

P6 = P 7 = 1 + sin(3A), (25) 

P8 = P 9 = sin A -I- cos(2A), 

77 = —2cos(4A). 

The dense (resp. dilute) phase is obtained for the pa¬ 
rameter range 0 < A < | (resp. f < A < |), and the 
corresponding CG coupling constant is p = ^ S (0, 2]. 

The loop model which is equivalent to the Q-state 
Potts model [20] corresponds to another choice of inte¬ 
grable weights: 


Pi = • • • = P7 = 0 , 

Ps = P9 = 1, (26) 

n = \/Q = —2cos(7rp), 

so that each edge is covered by a loop (complete packing). 
The CG coupling g S (0,1] covers only the dense phase 
in this case. 
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Transfer matrix method for measuring (7(01, 02 , da) 
General setup 

Let ($ 1 , $ 2 , ‘i’ 3 ) be three quasi-primary operators of a 
CFT, for which we want to compute the structure con¬ 
stant (7(<i>i, $ 2 ,‘ha)- The idea of our method is to con¬ 
sider the problem on the cylinder, and use the opera¬ 
tor/state correspondence with operators and states nor¬ 
malized so that ‘hj|0) = (where |0) is the ground 
state) to write 

C($i,$ 2,$3) = ($i|$2|$3)- (27) 

In order to achieve a finite-size estimation of this expres¬ 
sion, we simply need two ingredients: (i) the eigenstates 
|$i)l and |‘h 3 )L of the transfer matrix for the lattice 
model on the cylinder of circumference L, which con¬ 
verge to |$i) and |$ 3 ) as L —>■ 00 ; and (ii) the average 
value of the finite-size operator $ 2 ^^ between two basis 
states (a I and |6). The latter operator will typically scale 
like ~ < 1 ) 2 , where A /2 is a non-universal 

factor. This undetermined normalization is then easily 
eliminated by taking the ratio: 


within machine precision. This allows us to evaluate each 
scalar product on the left-hand side of (28) in just 0{Dl) 
operations, while using essentially the same amount of 
memory as in the direct method. 


Case of three purely electric operators 


Let us now explain how the above method is applied to 
the numerical computation of (7(di, 0 : 2 , ds). Our trans¬ 
fer matrix acts on the space of non-crossing link patterns 
between L points, where some arcs can pass through the 
periodic boundary condition horizontally (see Fig. § 1 ). 
The sector of the transfer matrix defined by giving a 
weight rii to non-contractible loops contains (the discrete 
analog of) the Verma modules of \Va.) with 


Oii G 


Q 

2 


Ci -|- 2Z 
26 


(29) 


where Ci = Arccos(ni/2)/7r. So we pick a state 
in the sector of weight m, for i = 1,2,3. The opera¬ 
tor associated to loop weight 712 , and sitting at 

site j G {1,..., L}, is defined through its average value 
between two basis states: 


^(<l)l|<l)^^^|<i)3)L 

L(‘h2|‘k^"'^|0)L 


L — 


(7($l,$2,«'3)- 


(28) 


The above relation gives a direct way to access the con¬ 
stant (7(<l)i, < 1 ) 2 , ‘f’ 3 ), by evaluating the eigenstates |<l)i) l, 
|< 1 ) 3 )l and |0 )l, acting with and performing the re¬ 
quired scalar products. The convergence in L is typically 
algebraic. 




where, after gluing the states a and b (see Fig. [5>), £{a, b) 
(resp. £ 2 ( 0 , b)) is the number of contractible loops not 
surrounding j (resp. surrounding j), and £i(a,b) (resp. 
£3(0,6)) is the number of non-contractible loops below 
j (resp. above j). The lattice operator contains 

contributions from the ’s with 0:2 defined by (291 with 


7 = 2, and their descendants under the Virasoro algebra: 


Time-efficient version 

In loop models, the most time-consuming part in this 
calculation is computing the scalar product, which takes 
0{D\) operations, where Dj^ is the dimension of the 
Hilbert space for the transfer matrix, growing exponen¬ 
tially in L. This is because the scalar product (a|6) be¬ 
tween two basis states is always non-zero. In the case 
when both |<I>i) and |<I> 3 ) are the lowest-energy states in 
two symmetry sectors of the transfer matrix, the method 
can be improved as follows. Denoting tj the transfer ma¬ 
trix in the sector where |<l)j)L lives, and \aj)L a basis 
state in this sector with L{aj\^j)L ^ 0, we use the fact 
that 


L{ai\tf^^2^t^\az)L 
\/L{ai\tV^\ai)L L(a3|ti^|a3)L 


M-f-oo 


L(‘hi|‘h^'^^|<k3)L 


for finite cylinder of 2M rows, where M ^ L. In practice, 
due to the exponential convergence with respect to M, 
only M ^ lOL iterations are needed to reach the limit 


E ■^(e 2 )T- 2 ^(^=)HQ/ 2 +e'/( 28 )+desc. 

e'o Ge2 +2Z 

(31) 

where A/'(e 2 ) is an undetermined prefactor. Thus, as in 
(28), the ratio 


converges to (7 (q;i, 02 , ds) for 62 G [—1,1]. The value of 
C for a charge 62 = 62 -I- 2k with k Gh can be obtained 
by picking the correct state \ Vafi)L and selecting the sub¬ 
dominant term with a prescribed scaling in 


the numerator and denominator of (28). 


Algorithm without scalar products 

The time-efficient version of our algorithm thus com¬ 
putes the partition functions Zijk = (d, 7 * 2 ,£* 3 ) 
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FIG. 5: (a) Link patterns for a transfer matrix of size L = 
4. We illstrate the Potts case (complete packing); the 0(n) 
model further admits empty points (dilution). (b) Gluing link 
patterns a and b is defined by juxtaposing a and the reflection 
of 6 in a horizontal mirror. The loop weights are indicated 
with the same colors as in Fig. 
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FIG. 6: Gylinder geometry with L = 4 and M = 3. The 
dashed lines are identified by the periodic boundary condition. 


on a finite cylinder of the square lattice, with circum¬ 
ference L and height 2 M rows, in the limit M ^ L. 
The three marked points reside on lattice faces, with Fi 
(resp. F3) situated just below (resp. above) the bottom 
row (resp. top) row, and F2 in the middle (see Fig. |^. 
For the Potts model we must take both L and M even to 
respect the sublattice parities in the mappings from the 
spin model [20] i whereas for the 0(n) model there are no 
such parity effects. 

The states of the transfer matrix are of the type shown 
in Fig. but with the possibility of dilution (empty 
sites) in the 0 (n) model. The boundary condition below 
the bottom row is taken as the first state in Fig.[^; in the 
Potts case this corresponds to free boundary conditions 
for the Potts spins. A similar boundary condition (hor¬ 
izontally reflected) is applied above the top row. Note 
that with these boundary conditions neither Fi nor F3 
can be be surrounded by a loop. The structure constant 
(^(di, 02, 0^3) is then computed from Zijk as explained 
around (101 in the main text. 


To build the lattice efficiently by the transfer matrix, 
we factorize the latter as a product of sparse matrices R, 
each corresponding to the addition of one vertex. The 
states are organized in a hash table. Each arc in the 
states (see Fig. carries two binary variables A^i2 = 0,1 
(resp. N23) that count the number of times modulo 2 that 
the arc has intersected the cut C12 (resp. C23) running 
from Fi to F2 (resp. from F2 to F3). The first stage in 
building a row with the transfer matrix is the insertion 
of the auxiliary space, i.e., the first horizontal edge that 
crosses C12 (resp. C23) on the lower (resp. upper) half 
cylinder. If that edge carries an arc, we initialize the 
corresponding A^i2 = 1 (resp. N23 = 1 ). In all other 
stages of the transfer process, when two arcs are joined 
their variables A^i2 and N23 add up modulo 2 . Finally, 
a loop resulting from the closure of an arc is attributed 
the weight nt given by ([^ in the main text (cf. Fig. 


Finite-size extrapolations 

In Figs. Hi and [^ we have displayed only the raw re¬ 
sults ( |I 0 | ) for sizes L = 4 , 5 ,..., 10 . The convergence to¬ 
wards the analytical result—which is already convincing 
by visual inspection—can be further improved by com¬ 
puting finite-size scaling (FSS) extrapolations 

C{L) = C{oo)+-fL-’^ ( 32 ) 

from three consecutive sizes L. For instance, in the case 
of Fig.j^the FSS exponent k was found to increase mono- 
tonically with rii —from 1.1 to 2.0 (resp. 0.6 to 2.1) in 
the dense (resp. dilute) case, for the parameter values 
shown—except for = 1 where (7 = 1 exactly for any L. 
The resulting agreement between ( 7 (oo) and ([^ amounts 
to at least four significant digits over the whole parameter 
range. 


The factor \/2 in the Potts model 

The correspondence between the lattice discretization 
and the continuum limit is somewhat singular in the 
case when non-contractible loops on the cylinder get zero 
weight (i.e., ni = 0 and/or 113 = 0 ). The states on which 
the transfer matrix acts are of two types, or parities, de¬ 
pending on whether the number of arcs crossing the pe¬ 
riodic boundary condition is even or odd. For instance, 
in Fig. the three states in the first (resp. second) line 
are even (resp. odd). The two subsets are related by a 
shift by one lattice spacing, and thus have the same car¬ 
dinality; this relationship corresponds to the Kramers- 
Wannier duality of the Potts model. 

It is easy to see [211 [2S] that when non-contractible 
loops are forbidden the parity just defined is conserved 
by the transfer matrix. In more mathematical terms. 
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the transfer matrix acts on a representation of the pe¬ 
riodic Temperley-Lieb algebra which is generically irre¬ 
ducible, but breaks down into two irreducible, isomorphic 
sub-representations when n* = 0. The lattice procedure 
we use restricts to just one of these sub-representations, 
while by continuation of the general Ui case, the the¬ 
oretical formula for the three-point function should in¬ 
volve, in the lattice regularization, the sum of both sub¬ 
representations. Associate formally to each of these rep¬ 
resentations two fields with the same conformal weight 


(but vanishing cross two-point function) Vi and V 2 . The 
theoretical formula, considering for instance the case 
where the three n* = 0, applies formally to 

(the ■\/2 factor arises from the normalization of the two- 
point function). In obvious short-hand notations one has 
VVV = jf transfer matrix calculation we 

v2 

measure ViViVi, we see that the result should be \/2 
larger than the theoretical value (which corresponds to 
VVV), in agreement with the numerical observations. 



